Unveiling First Order CMR Transitions in the Two-Orbital Model for Manganites 
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Large scale Monte Carlo simulation results for the two-orbital model for manganites, including 
Jalm- Teller lattice distortions, are here presented. At hole density x = 1/4 and in the vicinity of the 
region of competition between the ferromagnetic metallic and spin-charge-orbital ordered insulating 
phases, the colossal magnetoresistance (CMR) phenomenon is observed with a magnetoresistance 
ratio ~ 10, 000%. Our main result is that this CMR transition is found to be of first order in some 
portions of the phase diagram, in agreement with early results from neutron scattering, specific heat, 
and magnetization, thus solving a notorious discrepancy between experiments and previous theo- 
retical studies. The first-order characteristics of the transition survive, and are actually enhanced, 
when weak quenched disorder is introduced. 



Introduction. The manganese oxides known as man- 
ganites continue attracting considerable attention due to 
the presence of several competing tendencies in their rich 
phase diagrams, with a variety of spin, charge, and or- 
bital orders [l|, |2j. Moreover, these compounds display 
the famous CMR effect, where the resistivity is drasti- 
cally reduced by fields of a few Teslas. Early theoreti- 
cal investigations [3] remarked the importance of phase 
competition to understand the CMR, which occurs when 
an insulating state, typically also spin-charge-orbital or- 
dered, is close in energy to the low-energy ferromagnetic 
(FM) metallic ground state induced by double exchange. 
These effects were clear in simplified phcnomcnologi- 
cal models and resistor-network simulations when in the 
presence of disorder [4]. However, it is important to ver- 
ify at a more fundamental level if the basic model Hamil- 
tonians that are widely perceived as being realistic for 
manganites, including double-exchange, superexchange, 
and lattice distortion tendencies, are indeed compatible 
with the CMR phenomenology. Recently, our group and 
others initiated this effort based on Monte Carlo (MC) 
simulations [EHS)]. This is a challenging task since com- 
puter efforts based on the exact solution of the fermionic 
e g sector, for fixed MC-generated classical £ 2 g spin and 
lattice configurations, have a CPU time growing as N , 
with N the number of cluster sites times the number of 
orbitals. For this reason, recent efforts at realistic elec- 
tronic densities have used only one e g orbital, consider- 
ably reducing the CPU demands 8]. However, it is clear 
that the more complete two e g orbitals model must be 
investigated for a full understanding of the CMR physics. 
Alas, since this is far more time consuming, MC results 
for two orbitals are only scattered in the literature and 
they have not reached sufficient accuracy to unveil the 
true properties of this model in the CMR regime. 

In addition, an important qualitative difference still 
persists with regards to the order of the CMR tran- 
sitions. Until now, MC simulations carried out at 



CMR realistic hole densities, such as x — 0.33 or 
0.25, have indicated the presence of a rapid crossover, 
yet continuous, transition from the low-T FM metal 
to the high-T paramagnetic (PM) insulator. This oc- 
curs even in the clean limit of the MC simulations, 
i.e. without quenched disorder. While these results 
are in agreement with the phenomenology of many 
CMR manganites [2|, there are notorious examples 
where the CMR transition is of first order, such as for 
Pro.55(Cao,75Sr . 25 )o,4 5 Mn0 3 [§], Smn^SrnasMnOa [10, 
E3, (Ndo.5Smo.5)o.55(Bao.5Sro.5)o.45Mn0 3 [B], and oth- 
ers. In fact, early experimental studies already reported 
signs of irreversibility at the CMR transition for one of 
the most widely studied manganites, Lai-zCa^MnOs 
(LCMO) at x = 0.33, suggesting that its transi- 
tion is weakly first-order jl 31 ] - Subsequent magnetiza- 
tion and specific heat analysis clarified that indeed the 
x = 0.33 LCMO transition is of first order Q. Re- 
cent investigations have also addressed the order of the 
CMR transition, revealing multicritical characteristics in 
RE .5 5 Sro.45Mn03 (RE = Rare-Earth) [H|. 

Our main goal in this publication is to solve this puz- 
zling theory-experiment disagreement by carrying out a 
large-scale computational study of the two-orbital model 
for manganites at x = 1/4. It will be shown that our MC 
data unveils the presence of a first-order CMR FM-PM 
transition in this realistic model and density. 

Model and techniques. The two-orbital model Hamil- 
tonian used here is widely considered the minimal model 
for the proper description of the electronic properties of 
man gani tes and it has been extensively discussed be- 
fore [l6(. For this reason only a schematic descrip- 
tion will be here provided. The model is given by 
H = H DE + Hj T + H AF + H Dis 47]. At every Mn site 
it contains e g electrons and £ 2 g localized spins. -Side rep- 
resents the nearest-neighbors double exchange (DE) e g 
electron hopping at infinite Hund coupling. This term 
favors ferromagnetism away from the electronic density 



n = 1. Hjt is the coupling between fermions and oxygen 
breathing and Jahn- Teller lattice distortions, with cou- 
pling strength A. The sum H UE + H JT is Eq.(4) of Ref.0, 
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with the hopping along x for orbital a 
the energy unit. Haf = Jaf S(ij) ^i ■ Sj is the standard 
nearest-neighbors Heisenberg antifcrromagnetic coupling 
among the (classical) t2 g spins. Finally, i?Dis = J2i ^i^i 
is the quenched disorder term, with Ai = A x n (A = 
disorder strength; r\ — bimodal random number 1 or -1). 

Details of the computer simulation. The methodology 
of our MC simulations, based on the fermionic sector 
exact diagonalization described before, is standard and 
readers should consult |2j for details. However, the exten- 
sive characteristics of the present MC simulations merits 
a detailed discussion. The procedure was the following. 
For the clean-limit results A = 0, a random spin and 
lattice configuration was chosen at the highest studied 
temperature, T = 0.33, to initiate the runs. After 10,000 
MC steps, T was reduced and another 10,000 steps were 
performed. This "cooling down" process continued, us- 
ing a grid with 23 temperatures (lowest 1/300) that was 
denser near the critical temperatures To- After this al- 
ready demanding first step, the cooling-down results (ob- 
tained by measuring during the last 5,000 MC steps per 
T) revealed sharp, yet continuous, transitions at all cou- 
plings. However, at particular Ts in the MC time evo- 
lution indications of insufficient convergence were found. 
Thus, next for each T the results were further refined 
using additional 10,000 MC steps for extra thermaliza- 
tion, followed by 100,000 MC steps for measurements. 
By monitoring the results during the final MC evolution 
it was observed that this large effort produces now fairly 
stable results, revealing the first-order transitions in the 
region of phase competition discussed below. 

The clean-limit effort needed standard computer clus- 
ters with ~ 100 nodes. However, a similar procedure 
with quenched disorder, requiring 40 n disorder config- 
urations, 6 disorder strengths (A), and 20 temperatures 
for each Jaf, would have been impossible. Thus, the 
results with disorder reported below were obtained us- 
ing the UT-ORNL Kraken supercomputer (Cray XT5), 
where up to 5,000 processors were employed simultane- 
ously for periods of 24 hours. This amounts to a total 
computational effort in Figs. 3(c) and 3(d) (see below) of 
^250,000 hours (^30 years if ran serially). 

Clean limit phase diagram. The clean-limit phase di- 
agram of the two-orbital model is in Fig. [1] The chosen 
coupling A = 1.3 is representative of the regime where 
metallic and Jahn- Teller distorted insulating phases are 
in competition. The hole density x — 1/4 is in the real- 
istic range of experimental CMR investigations. Varying 
Jaf, and via the MC procedure described before, a sys- 
tematic analysis of spin and charge correlations (e.g. see 
Fig. [2]) leads to a phase diagram where the FM metallic 
phase generated by the DE mechanism competes with a 
previously discussed C0.25E0.75 insulator [18| (shown in 
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FIG. 1: (Color online) Clean-limit MC phase diagram of the 
x — 1/4 doped two-orbital model at A = 1.3 and on an 8 x 8 
lattice. Critical temperatures are estimated from spin struc- 
ture factors (Fig. 2) and real-space spin-spin correlations (not 
shown). First (second) order transitions are indicated with 
open (filled) symbols, and the shaded region is where the 
CMR is observed (Fig. 3). Also shown are the T—Q arrange- 
ments of classical spins and electronic charges (magnitude 
proportional to the radius of the circles) for the competing 
metallic and insulating states. Inset: total energy vs. Jaf. 
The straight lines crossing indicates a low-T first-order tran- 
sition. The spins were frozen to the FM and CE perfect spin 
states and the oxygen distortions were MC relaxed. 

Fig.QJ. This "CE" state is the natural generalization to 
x = 1/4 of the well-known x — 1/2 CE state llmj, and 
they only differ in the shape of the zigzag chains. Like 
in all computational based efforts, here relatively small 
clusters are used, but it is well known that locating the 
temperature range where the relevant correlation lengths 
are as large as the cluster size provides qualitatively cor- 
rect estimations of trends and critical temperatures. 

The main novelty of the phase diagram in Fig. [1] is 
the identification of a first-order transition separating the 
FM metallic phase from the PM high-T state, in the cou- 
pling range close to the competing CE insulator. Previ- 
ous MC investigations had not reached sufficient accu- 
racy to detect this first-order transition at realistic hole 
densities, such as x = 1/4. While previous efforts had 
clearly established the first-order nature of the direct low- 
T FM-CE transition (see inset of Fig. [IJ, now observing 
the more subtle first-order FM-PM transition represents 
qualitative progress in the modeling of manganitcs. 

Figure [^a) illustrates how the spin structure factor 
varies with T. At q=(0,0) and close to the region of CE 
competition, a discontinuous transition exists between 
the FM and PM states. There is also a first-order tran- 
sition from the CE to the PM state, close to the FM re- 
gion, our new results reveal that the true phase diagram 
of clean-limit manganite models actually has "multicrit- 
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FIG. 2: (Color online) Spin structure factor (5(q)) vs. T 
for the momenta corresponding to the (a) FM and (b) CE 
states. S(q) was used to determine the critical tempera- 
tures in Fig [T] For the CE phase, the relevant momenta 
are q=(ir/4, 3tt/4), (tv/2,tt/2), and (37r/4,ir/4), as well as 
q'=q+7r (i.e. rotating the zigzag chains by 90-degrees). 

ical" characteristics, with the robust FM-CE first-order 
transition at low T splitting into still first-order FM-PM 
and CE-PM transitions with increasing T, each ending at 
critical points (whose precise location is beyond our ac- 
curacy) . These results are compatible with recent multi- 
critical characteristics revealed in some manganites |15j . 

CMR effect. The first-order nature of the transitions 
observed here dramatically affects the transport proper- 
ties. Using the standard Kubo formula to calculate con- 
ductances [9, 0], Fig. [3] shows the resistivity (p) vs. T 
obtained in our simulations. Fig. [3Ja) are results in the 
clean limit, varying Jaf in the region of FM-CE com- 
petition. At couplings such as Jaf=0.148, p is insulat- 
ing upon cooling, closely following results for Jaf=0.16 
with a CE ground state. However, at the FM transition 
r c -0.02 (~ 100 K, if i aa =0.5 eV), p becomes metallic 
via an abrupt discontinuity, in excellent agreement with 
several experiments [9l4l4|. As Tq increases, moving the 
system further away from the CE state, the transition 
becomes continuous and at Jaf =0.12 the FM transition 
is barely noticeable in the slope of p vs. T. For complete- 
ness, in Fig. @Ic) results using a 12x12 cluster for only 
one set of couplings A- Jaf (due to its high CPU cost) are 
shown, indicating that cluster size effects are small [19j . 

Figure [3jb) contains the p curves in the presence of 
magnetic fields H. The observed trends are again in ex- 
cellent agreement with experiments |9l4l4 1 , with an over- 
all rapid decrease of p and with p-peak positions moving 
to higher T with increasing H, and all the curves merging 
at approximately room T (i.e. T ~ 0.06 if t aa ~ 0.5 eV). 
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FIG. 3: (Color online) (a) p vs. T, illustrating the metal- 
insulator transition varying Jap. Arrows indicate the mag- 
netic Tqs. Error bars are small (shown only at Jap = 0.148). 
(b) Magnetic field effects (along the z-direction) on the p vs. 
T curves, showing the CMR effect. (c,d) Quenched disorder 
influence (averages over 40 disorder configurations) on the p 
vs. T curves for (c) J A f=0.148 and (d) J A f=0.144, with first- 
order and continuous (but rapid) transitions, respectively. 

At a small field H = 7 x 10 _3 £ aa , the magnetoresistance 
(p(0) - p{H))/p{H) x 100% is - 10,000% also in good 
agreement with CMR phenomenology. 

Influence of Quenched Disorder. Experiments and the- 
oretical calculations have shown the importance that 
quenched disorder has over the CMR effect [2j. It is 
expected that the clean-limit fine tuning of couplings 
needed to obtain a CMR (Fig. [T]) will be removed once 
disorder is incorporated. Thus, it is important to ana- 
lyze the influence of disorder on our results. Using the 
on-site quenched disorder form described before, results 
are in Figs. [3fc,d). Panel (c) illustrates how the clean- 
limit first-order transition is eventually rendered contin- 
uous by increasing the disorder strength. However, the 
discontinuity in p first increases with increasing A before 
it is reduced. This result is compatible with the observed 
multicritical behavior even with disorder [15j . Also as in 
experiments, when the clean-limit transition is second or- 
der, quenched disorder decreases the T where the peak 
occurs while increasing p (panel (d)). 

The CMR effect observed in our clean-limit MC simu- 
lations can be understood qualitatively via measurements 
of a variety of observables, similarly as in previous inves- 
tigations [8[. For instance, Fig.|4ja) contains the spin and 
charge structure factors, S(q) and n(q), at the momenta 
of relevance for the CE phase, but in the CMR regime 
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FIG. 4: (Color online) (a) Spin S(q) and charge n(q) struc- 
ture factors, shown together with the resistivity p and the 
Fermi-level inverse density-of-states, vs. T. The "q=CO" 
shorthand stands for the characteristic momenta of the CE 
state (q=(-7r/2,7r/2), (n,n), and (37r/2, 37r/2), and momenta 
obtaining by rotating the zigzag chains by 7r/2). n(q) is ob- 
tained by Fourier-tranforming the quantum correlations be- 
tween local charge densities among e g orbitals x 2 — y 2 and 
3z 2 — r 2 , and averaging over the MC steps, (b) Pseudogap in 
the density of states N(uj) as T is reduced, (c) p vs. T for a 
12x12 cluster, showing also a first-order transition. 

where the T = ground state is FM. As T decreases, not 
only p rapidly increases but so do S(<ice) and n(qco), 
showing that the system behaves as if the ground state 
were CE, developing robust CE short-range correlations. 
However, at Tc an abrupt transition occurs to the true 
FM ground state. This switch from CE-dominated to 
FM-dominated characteristics with cooling may occur if 
the high-T short-range-ordered CE state has a high en- 
tropy. Note also that at low-T, the charge correlations in 
the FM state are still robust at short distances. Finally, 
also note that at a T*/t aa ~ 0.07 — 0.08, considerably 
higher than Tc, the CE tendencies start developing upon 
cooling [20] , and this occurs concomitantly with the pres- 
ence of a density-of-states pseudogap (Fig. Bib)), as also 
observed in photoemission experiments |21| . It is grat- 
ifying to observe similar results above Tc for both the 
present model and the model studied in Ref. 8. Finally, 
note that the many observations of FM signals above Tc 
in several manganites [2J are not incompatible with the 
clean-limit Fig. U since quenched disorder is known to 
increase the strength of the FM component. 

Conclusions. The observation of first-order CMR tran- 
sitions in models for manganites was here reported, solv- 



ing a notorious theory-experiment discrepancy. A large- 
scale computer simulation was needed to reach our con- 
clusions. Robust CMR ratios were found, as well as a 
CMR state above Tc with short range CE characteris- 
tics. Weak quenched disorder preserves the first-order 
transitions, as in recent experiments 15j. 
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